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We reinvestigate the Deterministic Lattice Gas introduced as a paradigmatic model of the 1/f 
spectra (Phys. Rev. Lett. 64, 3103 (1990)) arising according to the Self-Organized Criticality 
scenario. We demonstrate that the density fluctuations exhibit an unexpected dependence on 
systems size and relate the finding to effective Langevin equations. The low density behavior is 
controlled by the critical properties of the gas at the absorbing state phase transition. We also 
show that the Deterministic Lattice Gas is in the Manna universality class of absorbing state phase 
transitions. This is in contrast to expectations in the literature which suggested that the entirely 
deterministic nature of the dynamics would put the model in a different universality class. To our 
knowledge this is the first fully deterministic member of the Manna universality class. 
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I. INTRODUCTION 

The aim of statistical mechanics is to establish a de- 
scriptive bridge between the micro world and the emer- 
gent observables at the macroscopic level. To com- 
prehend the macroscopic world we make use of a few 
collective variables such as density or current. Ideally 
equations of motion for these variables should be rig- 
orously derived from the equations of motion of the 
components, however this is seldom possible even for 
the simplest systems. An effective approach consists of 
semi-phenomenological Langevin equations It is well 
known that although Langevin's approach is of great in- 
tuitive appeal it suffers from the fact that the details of 
the stochastic driving are very difficult to derive. Here we 
consider a microscopically deterministic model, namely 
the Deterministic Lattice Gas (DLG)0] and show how 
a bulk noise term, absent for small system sizes 0, 
is generated as the system size increases. In terms of 
driven diffusion equations it means that the behavior of 
the power spectrum of the density fluctuations can be 
understood in small systems as arising from a bound- 
ary driven diffusion equation^, Q without any Langevin 
bulk noise term. However as the system size is increased 
it turns out that a bulk noise term is needed in order to 
generate the correct power spectrum. The situation is 
particularly interesting because the model exhibits two 
regimes of very different behavior. At large density the 
fluctuations are well described by the appropriated diffu- 
sion equation. At low density the nearby absorbing state 
phase transition controls the fluctuations, which turns 
out to be in the Manna universality class. 



The DLG is a model of a gas of particles undergo- 
ing overdamped dynamics. Numerical simulations and 
analytic studies concluded that the fluctuations in the 
number of particles exhibit robust 1/f fluctuations in 
several configurations. In the context of Self-Organized 
Criticality (S0C)|1 this is of considerable interest since 
the original motivation for SOC was to establish a simple 
paradigm that might explain the widespread observation 
of 1// temporal fluctuations and spatial fractals. DLG 
is thus an example of this scenario in contrast to the 
sandpile model used by Bak, Tang and Wicscnfeld to in- 
augurated which turned out to contain no 1// spectraQ . 



II. MODEL 

In the DLG particles interact through a nearest- 
neighbor repulsive central unit force and double occu- 
pancy is not permitted. All particles are updated si- 
multaneously by moving the particles deterministically to 
neighbor sites according to the vector sum of the forces 
they are subject to. If two particles want to move to 
the same site, the particle subject to the strongest force 
is moved, while in the case of equal forces no particle is 
displaced. Periodic boundary conditions are considered 
and we consider the model in dimension d = 2. 

The number of particles N{t) in a sub- volume of the 
lattice exhibits interesting temporal fluctuations. We de- 
termine the power spectrum S{ f) oi N{t) from the square 
of the absolute value of the Fourier transform. Successive 
time sequences are averaged in order to achieve sufficient 
statistics: 
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The angular brackets denote averaging over many differ- 
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ent time series. 

The spectrum S{f) was shownQ to satisfy S{f) ~ 
1//^ at the particle density p = 0.3, while previous works 
showed that the same result is obtained in a wide range of 
densities with a drive at the boundary 0. In those papers 
the maximum linear system size considered is L = 128, 
we will see that different behaviors are observed for larger 
systems. In Fig. [T]we show how the effective exponent 
fi changes from fi — 1 to fi~S/2 as L increases. It is 
interesting to relate this change in the temporal fluctua- 
tions to the appropriate macroscopic equation of motion 
for the particle density p{r,t). 
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FIG. 1: Scaling behavior of the spectrum <S'(/) oc /^'' of the 
total number of particles N{t) in the DLG for increasing linear 
sizes of the lattice L. A crossover from /i = 1 for small L to 
/i = 1.5 for large L is observed. Particle density p = 0.5. 



Simulations show that although the particles move de- 
terministically in the DLG model they do behave like ran- 
dom walkers in the sense that their square displacement 
is hnear in time. This suggests that p{r,t) evolves ac- 
cording to a diffusion equation and integrating the diffu- 
sion equation allows us to extract the temporal evolution 
of N{t) = Jy drp{r,t), where V denotes the measuring 
sub-volume. Assume that p{r, t) satisfies the diffusion 



dp 
dt 



(2) 



It has been shownj^Q that the power spectrum of N{t) 
depends on the details of the diffusion equation. If the 
equation is driven by white noise on the boundary of 
V and the noise term ^ is absent the power spectrum 
exponent derived from Eq. ^ is p = 1. If in contrast 
a conservative bulk noise term is included in Eq. ^ 
one obtains instead p = 3/2. Including non-linearities 
in Eq. ^ will not influence the behavior of p, see [1]. 
It is this observation that allows us to conclude that a 
conservative noise term is generated in the macroscopic 
Langevin description when the sub-volume V becomes 
sufficiently large. 



III. S{f) SCALING FOR LARGE LATTICES 

We now consider the DLG for much larger lattices than 
previously studied and our aim is to understand the fluc- 
tuations spectrum as a function of density. We plot S{f) 
in Figure [2] and find that the value of p depends on the 
density: 




p > 0.3 

p ~ = 0.245 



(3) 



where pc is the critical density of the absorbing phase 
transition (APT in the following). The observed expo- 



< 

CO 




0.01 



1e-05 



0.0001 



0.001 0.01 
f 



0.1 



FIG. 2: Scaling behavior of the spectrum S{f) oc of the 
total number of particles N{t) in the DLG for different par- 
ticle densities p. A crossover from /i ~ 1.8 at p ~ pc = 0.245 
to p = 1.5 at p ^ Pc is observed. S{f) has been multiphed 
by different constants for different densities p to visualize the 
scaling exponents properly. Lattice hnear size L = 1000. 



nent at high densities is the same of a gas of random 
walkers while at low densities the power spectrum 
scaling exponent is determined by the critical properties 
of the DLG at the APT, as we shall discuss in the fol- 
lowing. 



IV. ABSORBING PHASE TRANSITION 

At very low densities the DLG enters configurations 
in which all particles are far away from each other and, 
due to the short range interactions, the particles become 
unable to move, i.e. the dynamics is frozen. This is 
called an absorbing state phase transition between an 
active phase, above the critical density pc, and an inactive 
below. 

The control parameter of the transition is the particle 
density p, while the order parameter is the density of 
active particles pa (this is the average fraction of particles 
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that move in one time step). Near the transition: 

MM -If (4) 



with 5p = p — pc- The determination of the critical 
density pc is obtained by varying its value until one 
observes a straight line in a log-log plot, which gives 
Pc = 0.24500(2). The brackets indicate the statistical 
uncertainty on the last digit, which is estimated by look- 
ing at deviations from a straight line in the same plot. A 
regression analysis yields the value of the order parame- 
ter exponent /3 = 0.634(2). 
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FIG. 3: Density of active sites pa as a function oi Sp = p — 
Pc- The determination of the critical density pc is obtained 
varying its value until data points are aligned on a straight 
line in a log- log plot. Lattice linear size L — 1000. Error bars 
are smaller than symbols. 
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FIG. 4: Determination of the exponent 7'. Lattice linear size 
L = 1000. 



spontaneous creation of active particles, which prevents 
the system to enter an absorbing state and thus allows to 
study scaling properties at the critical point p ^ pc- Here 
we make use of the implementation of the external field 
for APTs that was developed in [l^. Since the control 
parameter p is a conserved quantity, the external field 
must preserve the conservation: at each time step, after 
performing the DLG update, we choose randomly hL^ 
particles on the lattice and move each of them to one of 
its empty neighbors. In this way inactive particles may 
be activated and the number of active sites is increased, 
while the density of particles p is conserved. The intro- 
duction of the external field prevents the system from 
falling into an absorbing state and the order parameter 
Pa docs not vanish at p = pc, as can be seen in Fig. [5] 



The fluctuations of the order parameter Apa = 
[{pD ~ (Pa)^] are plotted in Fig. |4l Approaching the 
transition point Apa increases and diverges at pc. Close 
to the critical point fluctuations scale as 



Apa ^ Sp 



(5) 



in the active phase. Using a regression analysis one gets 
7' = 0.40(2). 



A. External field 

A feature of finite size effects is that the system may 
pass within one simulation from one phase to the other. 
This is caused by critical fluctuations, that diverge while 
approaching the critical point: as soon as the correlation 
length is of the order of L the system may pass to the 
absorbing phase even in a small region above the critical 
point Pc- Similar to equilibrium phase transitions it is 
possible to apply an external field h conjugated to the 
order parameter, i.e., in our case the field corresponds to 
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FIG. 5: Density of active sites pa vs particle density p for 
different values of the external field h. At /i > the system 
does not fall into an absorbing state. Error bars are smaller 
than symbols. Lattice linear size L=1000. 



The order parameter and its fluctuations are assumed 
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to be generalized homogeneous functions of the control 
parameter p and of the applied field h: 

Pa{5p.h)^\-''R{5p\hy) (6) 

Apa{5p,h) = X^' b{5p\h\'') (7) 

where A > 0. From equation ([5]) one can obtain an esti- 
mate of the exponent a plotting pa against h in a log-log 
plot: 

Pa{5p^Q,h)^hP'- (8) 

This gives a — 2.19(1)^. An independent estimate of 7' 
can be obtained through equation ([7]) in the same way: 

Ap,(5p==0,/i)^/i-T> (9) 

which gives 7' = 0.365(7). We consider the weighted av- 
erage of the two determinations of 7' as our best estimate: 
7' = 0.37(1). 

Data collapses can be produced choosing A = h^^^'^ in 
equations © and ([7]), see figures IH and [71 

h-l'/'^PaiSp, h) = R{Sph-^/% 1) (10) 




FIG. 7: Data collapse as in equation In the inset data 

before the collapse are shown. Lattice linear size L = 1000. 



B. Finite size scaling 

Similar to equilibrium critical phenomena we assume 
that the system size L enters the scaling forms ([6|) and 
([7]) as an additional scaling field, i.e., 

PaiSp,h,L) - \-''Rpbci5pX,hy,L\-''^) (12) 



h'^'/'^ApaiSp, h) = biSph-^'", 1) 



(11) 



Apa{Sp,h,L) = AT DpbciSpX,hX'',LX-''^) (13) 
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FIG. 6: Data collapse as in equation (|10p . In the inset data 
before the collapse are shown. Lattice linear size L = 1000. 



The error on /3/(t is computed as the dispersion on its value when 
measured at pc,l = 0.24498 and pc,2 = 0.24502 (respectively the 
lower and upper bound on our determination of the critical den- 
sity pc = 0.24500(2)). Errors on other quantities are estimated 
in the same way, when possible. 



where the exponent i/± describes the divergence of the 
spatial correlation length, i.e., oc 5p~^^. The uni- 
versal scaling functions depend on the particular choice 
of the boundary conditions, although in the thermody- 
namic limit Rpbc{x,y, 00) ~ R{x,y) and Dpbc{x,y, go) = 
Dix,y). 



Following 11| we consider the fourth order cumulant 
Q, which is defined as: 



(14) 



For non-vanishing order-parameter the cumulant tends 
to Q = 2/3 in the thermodynamic limit. One expects it 
to obey the scaling form: 

Q{Sp,h,L) ^ QpU^pX^hX^LX-"^) (15) 

Choosing A = L^/'^^ in equation (|T5|) ai 5p ~Q we obtain 
the following equation: 

Q(0,/i,i) = Qpbe(0,/ii"/''^,l) (16) 

which enables us to determine v±_ through a data collapse 
by plotting Q{0,h,l) against hL"^^"^ as in Fig. [5] Best 
results are obtained for i'± = 0.83(5)^. 



^ The error is estimated looking at data collapses for different val- 
ues of the exponent. 
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FIG. 8: (lin-log) Data collapse as in equation (|16|) for the FIG. 9: Data collapse as in equation (|19p to estimate z. In 

determination of vx_ . In the inset data before the collapse are the inset data before the collapse are plotted. 

shown. 



C. Dynamical scaling 

Starting simulations of the DLG from a random distri- 
bution of particles above the critical point pc the density 
of active sites pa decreases in time and tends to its steady 
state value. At the critical point p = pc the order param- 
eter decays algebraically as: 

Pa{5p^Q,h^Q,t)^t-'' (17) 

Simulating a lattice of linear size L = 4000 at p = pc we 
obtain a = 0.41(1). 

As can be seen in the inset of Fig. ^ a finite system 
size limits this power law behavior and one expects the 
following ansatz to hold at criticality: 

Pa{L,t) = X-"''iiR'^,,{tX-''ii,LX-^^) (18) 

Setting aijLX^^"^ = 1 one obtains: 

Pa{L,t)=L-'"%,,{tL-^-,l) (19) 

with z — v\\lv}_. Through equation (|19p one can ob- 
tain an estimate of z performing a data collapse, plotting 
L°''^ Pa against tL~'^ as in Fig. Inland varying the value of 
the exponent until the collapse is satisfactory. The best 
data collapse is obtained for z = 1.5(1)'^. 

D. DLG universality class 

After having determined the scaling exponents we are 
able to address the question concerning which universal- 
ity class the DLG belongs to. In table U we compare 



the measured values of the DLG critical exponents with 
those of the Manna universality class [12], showing that 
they are compatible with each other^. Some of the errors 
in our determinations of the exponent are significantly 
larger than the errors on the estimates of the Manna 
universality class available in the literature: this is prob- 
ably due to the fact that we cannot afford to perform 
as many iterations of the model (to compute ensemble 
average) as those performed for the CLG, CTTP and 
Manna model[l2|, since the DLG is computationally very 
demanding to simulate. It is also likely that stochas- 
tic models are more self-averaging than our deterministic 
model is and therefore allows more accurate estimates of 
the exponents in less iterations. The shape of the univer- 
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Manna 


0.639(9) 


0.799(14) 


1.225(29) 


2.229(32) 


DLG 


0.634(2) 


0.83(5) 


1.2(1) 


2.19(1) 




7' 
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Q 


z 


Manna 


0.367(19) 


1.590(33) 


0.419(15) 


1.533(24) 


DLG 


0.37(1) 


1.54(1) 


0.41(1) 


1.5(1) 



TABLE I: The measured critical exponents for the DLG and 
the correspondin g cr itical exponents for the Manna universal- 
ity class in d = 2p^. 



sal scaling functions is remarkably similar to those that 
can be found in the literature for the Manna universality 
class. This is seen e.g. by comparing our Fig. [7] with Fig. 
4 in [l3l, see also [H]. From the scaling exponents and 
the scaling functions we find compelling evidence that 
the DLG belongs to the Manna universality class. 



^ The error is estimated looking at data collapses for difEcrcnt val- 
ues of the exponent. 



* The estimate of is obtained from z = v\\lvi_, knowing z and 
V ^. The value of 7 is obtained through Widom law (7 = cr — 
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V. POWER SPECTRUM AT CRITICALITY 

In this section wc want to relate the exponent of the 
power spectrum at low density to the scaling properties 
of the DLG at the critical point pc of the APT[13|. As we 
will see the scaling behavior of the correlation function 
Ca(r,t) = (pa(r,t)pa(0,0)) - (pa (0,0)2) ^j^g critical 
density p — determines the power law exponent in 
the scaling of the spectrum Sa{f) of the total number of 
active particles Na{t) in the lattice. 
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In stationary directed percolation processes above crit- 
icality it has been observed [l^] that the correlation func- 
tion Ca(r, t) a.t t = first decays in space algebraically as 
j.-P/i^±^ until it saturates at a constant value at r > 
In the saturated regime the two sites become uncorre- 
lated so that this value is just equal to the squared sta- 
tionary density of active sites p^ . 

One can use this fact to extract the scaling exponent 
for Sa{,f), making use of the Wiener-Khinchin theorem: 

Ca{t) = {Na{t)Na{Q)) - {Na{Of) 

' [ d^r [ r' Pair, t) Pair', 0)) + ... 



V l^d'r pair,t)paiO,0)) + 



d^r Ca{r,t) + 



^1/^(2-/3/1/^) 



dr r r-f^/'^^Cait/r') 



so that: 



Saif) 



1 

2^ 



dt Cait) e 



f 



r 



(21) 



With our estimates of the 



with p = l + \[2-^ 

critical exponents for the DLG we find ji = 1.82(6), while 
with the best estimates of the Manna universality class 
exponents [T^ one finds p = 1.78(2). 

At low densities we expect that fluctuations in the total 
number of particles Nit) in the lattice are triggered by ac- 
tive particles, therefore we expect the two spectra £'(/) = 
^Jdt Nit) e-*2^/* and Saif) = ^ ! dt Nait) e-^^^f' 
to show the same scaling behavior at p ~ pc, as it is 
confirmed by simulations (see Fig. [T0| . 

Computing the power spectrum of the total number 
of particles in the box at p = 0.24506 ~ pc we observe 
a power spectrum compatible with the predictions, con- 
firming that fluctuations in the total number of particles 
in the sub-volume V are determined by the scaling be- 
havior of the correlation function at criticality. 



FIG. 10: DLG: Scaling behavior of the spectrum S{f) oc 
of the total number of particles in the box Nit) and of the 
spectrum Saif) of the total number of active particles near 
the critical density p = 0.24506 ~ Pc- The spectra have been 
multiplied by arbitrary factors to visualize the scaling expo- 
nent properly. Lattice linear size L — 1000. 



VI. CONCLUSIONS 

We have investigated the scaling properties of the 
power spectrum Sif) of the total number of particles 
in a sub-volume of the lattice for the DLG: our simu- 
lations show that the DLG is not characterized by 1/f 
fluctuations as had been previously observed, but present 
a much more variegated picture. At high densities the 
power spectrum is the same as observed in a gas of ran- 
dom walkers, while at low densities the spectrum scales 
as Sif) ^ 1//^'* and we have shown that this exponent 
is determined by the decay of the density-density corre- 
lation function near criticality. 

Furthermore we have shown that the deterministic lat- 
tice gas is in the same universality class as the Manna 
model[l3|- To our knowledge this is the first example 
of a completely deterministic and non-chaotic system in 
this universality class. Perhaps this is not too surprising 
given that we found that a stochastic bulk noise is gener- 
ated in the effective Langevin description as the system 
size is increased. From this point of view the DLG be- 
comes similar to the other members of the Manna class, 
which all involved a stochastic element in their micro- 
scopic dynamics. Examples of models in the Manna class 
includes the Oslo model [Tsl and the Manna model[l3|- In 
the Oslo model the local threshold for relaxation is up- 
dated stochastically at every relaxation and in the Manna 
model particles move to stochastically chosen neighbor 
sites. 

Further studies could focus on the critical behavior of 
the DLG in one, three and higher dimensions, or on deter- 
ministic variants of the Manna and Conserved Threshold 
Transfer models [l6j. to address the question of how a de- 
terministic micro-dynamics affects the critical behavior 
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of these systems. Recently a method |17| has been de- 
vised to check wether a model belongs to the Manna or 
to the DP universality class, that have similar critical ex- 
ponents. This method, that involves the study of critical 
properties near an absorbing or reflecting wall, could be 
used to confirm the fact that the DLG belongs to the 
Manna class. 



so that we expect the following ansatz to be satisfied at 
h = near the critical point pc'. 

DtiSp,h^O)^Sp^ (25) 

This is verified in Fig. Illl and a linear regression analysis 
gives /3 = 0.637(4). 



Appendix 
Non-universal metric factors 

When we express the order parameter and other quan- 
tities as generalized homogeneous functions like in equa- 
tions ([6]) or ([7]) we can introduce the non-universal metric 
factors Op, a a and so that once they are chosen in a 
specified way, the universal scaling functions are the same 
for all systems within the same universality class, while 
these factors arc non-universal in the sense that they can 
be different for models in the same class (see [l^l for fur- 
ther details): 

PaiSp, h) = \-PR{ap5p\, anhX") (22) 

aAAp,(,Sp, h) = X^'biapSpX, a^hX") (23) 

In table HI] we report the measured values of the non- 
universal metric factors for the DLG. 
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DLG 


1.741(9) 


27(2) 


0.028(2) 


40.1(1) 


CLG 


0.509 


9.241 


0.062 


11.22 


CTTP 


0.341 


45.42 


0.013 


24.90 


Manna 


0.211 


78.56 


0.007 


35.53 



TABLE IL The measured non-universal metric factors for the 
DLG and for the CLG, CTTP and Manna models^ (the 
uncertainty on these values is less than 5%), in 2d. 



Tracer diffusion constant as an order parameter 

It has been recently observed [l^ that in some one- 
dimensional stochastic sand-piles the tracer diffusion co- 
efficient Dt ((i?^(i)) = 2 Dtt, where R{t) is the displace- 
ment of a particle at time t) scales in the same manner 
as the density of active particles, so that it represents 
an alternative definition of an order parameter. Here we 
show that the same behavior is found in the DLG, which 
enables us to obtain independent estimates of the scaling 
exponents /? and a. 

Scale invariance implies that the new order parameter 
Dt can be written as a generalized homogeneous function 
(for an infinite system): 
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ln(p - Pc) 

FIG. 11: Determination of the order parameter exponent 
/3 through the scaling of the tracer diffusion coefficient: 
Dt{5p,h = 0) ~ ^p". L = 1000. 



Also, for zero field h at criticality one has (see Fig. [T2|) : 
Dt{Sp^O,h) hf^/" (26) 

and a linear regression analysis gives /3/ct = 0.302(4), 
therefore cr = 2.11(4). Both two exponents are compati- 
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FIG. 12: Determination of jS/a through the scaling of the 
tracer diffusion coefficient in the presence of an external field 
h: Dt{Sp = 0,/i) ~ /i"/'". Error bars are smaller than sym- 
bols. L = 1000. 



Dt{Sp, h) = X-'^f{SpX,hX'') 



(24) 



blc with the previously determined ones for the DLG. 
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